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1.  Introduction 


In  tliis  paper  we  focus  on  an  inverse  eigenvalue  problem  for  unitary  Hessenberg  matrices  with  pos¬ 
itive  subdiagonal  elements.  Throughout  this  paper  all  Hessenberg  matrices  are  upper  Hessenberg 
matrices.  This  class  of  matrices  bears  many  similarities  with  the  class  of  Jacobi  matrices,  i.e.  real 
symmetric  tridiagonal  matrices  with  positive  subdiagonal  elements.  Any  matrix  in  either  class  is 
normal  and  has  distinct  eigenvalues.  Both  n  X  ra  Jacobi  matrices  and  n  x  n  unitary  Hessenberg 
matrices  with  positive  subdiagonal  elements  can  be  parameterized  by  2ti  —  1  real  parameters.  This 
is  obvious  for  Jacobi  matrices;  for  unitary  Hessenberg  matrices  this  parameterization  is  described 
below.  Since  unitary  Hessenberg  matrices  with  positive  subdiagonal  elements  are  determined  by 
0(n)  parameters,  one  can  develop  efficient  algorithms  for  this  class  of  matrices.  These  algorithms 
are  analogous  with  algorithms  for  Jacobi  matrices.  For  example,  the  unitary  QR  algorithm,  de¬ 
scribed  in  [Gr2],  has  many  similarities  with  the  QR  algorithm  for  Jacobi  matrices.  Another  example 
is  provided  by  the  divide-and-conquer  methods  that  have  been  developed  for  both  the  tridiagonal 
and  unitary  eigenproblems  [Cu],  [DS],  [GRl],  [GR2]. 

In  this  paper  we  show  another  analogy  of  unitary  Hessenberg  matrices  with  Jacobi  matrices,  namely, 
that  a  unitary  Hessenberg  matrix  H  with  positive  subdiagonal  elements  is  uniquely  determined 
by  its  eigenvalues  and  the  eigenvalues  of  a  unitary  rank-one  perturbation  of  H .  The  matrix  H 
can  be  constructed  using  O(n^)  arithmetic  operations  using  an  inverse  unitary  QR  algorithm. 
Similar  results  for  Jacobi  matrices  are  well-established  [BGl],  [BG2],  [GH].  The  inverse  unitary  QR 
algorithm  is  analogous  with  the  algorithm  described  in  [GH]  for  Jacobi  matrices. 

2.  Unitary  Hessenberg  Matrices  and  Szego  Polynomials 

We  refer  to  a  finite  Schur  parameter  sequence  of  length  n  as  a  sequence  of  complex  numbers 
with  \'ij\  <  1  for  1  <  ^'  <  n  and  |7„|  =  1.  Also  define  the  complementary  Schur  parameters 
by  <Tj  :=  \/l  -  ITjI*-  Associated  with  the  finite  Schur  parameter  sequence  is  a  unitary 

Hessenberg  matrix  H  with  positive  subdiagonal  elements 

ft  =  •••  1  ^n-li  7n)  :=  Gi (71)^2(72)  .  •  .G„_i(7„_i)G'„(7n),  (2.1) 

where  the  Givens  reflector  is  the  identity  matrix  of  appropriate  size  except  for  the  2x2 

principal  submatrix 

with  the  bar  denoting  complex  conjugation.  The  matrix  on  the  right  in  the  product  (2.1)  is  defined 
by  Gniln)  ■=  diag[l,  1, . . . ,  1, -7n].  The  nonzero  entries  oi  H  =  [Vj,k]],k=i  ^ben  given  by 
Pj+ij  :=  (Tj  and  ijj^k  ■=  -Jj-iCTjcrj+i  ...ak-iJk  for  I  <  j  <  k,  where  70  :=  1. 

It  is  easy  to  see  that  every  n  X  n  unitary  Hessenberg  matrix  H  =  [’Ij, >  0  is 
uniquely  determined  by  a  finite  Schur  parameter  sequence  of  length  n.  In  fact,  the  Schur  parameters 
{7j}j_i  and  the  complementary  Schur  parameters  can  be  determined  from  H  by 

I  <  j  <  n;  !<><«• 

Hence,  we  have  a  one-to-one  correspondence  between  Schur  parameter  sequences  of  length  n  and 
nxn  unitary  Hessenberg  matrices  with  positive  subdiagonal  elements.  This  Schur  parameterization 
of  unitary  Hessenberg  matrices  with  positive  subdiagonal  elements  show’s  that  these  matrices  are 
determined  by  2n  —  1  real  parameters:  the  real  and  imaginary  parts  of  7j  for  1  <  J  <  n,  and  the 


j  i  +  1  _  -7j  (Tj 

j  J  +  f  J  [^3  J  ’ 


argument  of  ")„■  To  avoid  numerical  instability,  however,  we  also  retain  the  complementary  Schur 
parameters. 

Let  II  =  Hn  :=  //(■)!, ...  ,7n-i,7n),  and  let  Hk  :=  G’i(7i)G2(72)  •  •  ■<T/.-i(7fc-i)<i’fc(7/;)  be  the 
leading  principal  submatrix  of  II  =  Hn  of  order  k.  Introduce  the  functions 

0fc(A)  :=  eJiXl- Hk)-''eu 
Vi(A):=det(A/-^*), 

trfc(A)  :=det(A/-if;'_i), 

where  Ci  :=  [1, 0, . . . ,  0]^,  and  denotes  the  trailing  principal  submatrix  of  Hk  of  order  k  —  1. 

Thus,  d>/t(A)  =  7r*:(A)/^’*^( A).  The  following  proposition  can  be  verified  by  induction. 

Proposition  2.1.  The  polynomials  and  7rjt(A),  fc  >  0,  satisfy 


li'fc(A)  7rfc(A) 

X  Ik 

’V’fc-i(A) 

7rjt-i(A) 

V’o(A) 

1  1/A 

t\.(A)  fk(X) 

7k>^  1 

L  j 

^\-i(A) 

Ilk-iW 

) 

V-o(A) 

1  0 

It  follows  that  for  each  k,  the  polynomials  ipk{^)  ■—  are  obtained  by  reversing  and 

conjugating  the  coefficients  of  Vjt(A).  From  the  initial  conditions  we  recognize  t/’fc(A)  to  be  the  I’th 
S:ego  polynomial  determined  by  the  Schur  parameter  sequence  □ 

The  Szego  polynomials  {'ibk}'k=o  orthogonal  with  respect  to  a  discrete  measure  on  the  unit 
circle.  This  measure  assigns  a  positive  weight  Uk  to  each  zero  \k  of  These  weights  are  the 

numerators  in  the  partial  fraction  decomposition 


it=i 

see  [Grl]  for  details. 

Let  H  =  UXU",  with  A  =  diag[Ai,  A2, . . . ,  A„]  and  U  unitary,  be  the  spectral  resolution  of 
H  =  /f  (71 , . . . ,  7„_i ,  7„  j,  where  +  denotes  transposition  and  complex  conjugation.  Let  u  :=  U^ei 
be  the  vector  containing  the  first  components  of  the  eigenvectors  of  H .  Since  H  has  nonzero  sub¬ 
diagonal  elements,  every  entry  of  u  =  [ui,U2»  •  •  Vn]^  is  nonzero,  and  we  normalize  U  so  that  each 
t’fc  >  0.  Then 

n  o 

<?i„(A)  =  u*(A/-  A)-»u  =  Yi  T^- 

k=l 

Hence,  the  weights  Uk  =  v\  are  guaranteed  to  be  positive,  and  Yl'k=.\'^k  =  !• 

3.  The  Inverse  Unitary  QR  Algorithm 

Given  n  distinct  unimodular  complex  numbers  {A;t}Jt=i  associated  positive  weights 
we  can  construct  a  unitary  Hessenberg  matrix  H  with  the  Xk  and  Vk  equal  to  the  eigenvalues  and 
first  components  of  the  corresponding  eigenvectors  of  //,  respectively.  This  construction  is  achieved 
using  an  inverse  QR  algorithm,  which  is  analogous  with  the  procedure  of  [GH]  for  real  symmetric 
tridiagonal  matrices. 

The  required  Hessenberg  matrix  is  obtained  by  performing  a  sequence  of  elementary  unitary  simi¬ 


larity  transformations  to  transform  the  matrix 


u 

A 


to  a  Hessenberg  matrix 


S 


^1 

H 


without 


using  or  changing  the  arbitrarj'  entry  6.  Then  H  =  7n-i  >  7n)  has  the  desired  eigenvalues 

and  associated  eigenvectors. 

The  idea  is  to  build  the  Hessenberg  matrix  by  adding  weight-abscissa  pairs  one  at  a  time.  Suppose 
that  we  have  constructed  the  unitary  Hessenberg  matrix  Hm  •=  .^'(71.  •  •  •  >7Tn-ii7m),  for  some 
m  <  n,  corresponding  with  the  weight-abscissa  pairs  .  Let  oq  :=  (53^=1 

assume  that  the  first  components  of  the  eigenvectors  of  Hm  are  {(t»*/(7o)}”=i  •  In  order  to  add  the 
weight-abscissa  pair  (f*,  A)  and  construct  the  corresponding  (m  -f- 1)  X  (m  -I-  1)  unitary  Hessenberg 
matrix  :=  //(AJ,  A.^, . . .,  we  perform  a  sequence  of  unitary  similarity  transformations 

to  put  the  (m  +  2)  X  (w  +  2)  matrix 

'  6  V  ao 
V  A 

8  V  <7oel]  (To  X  X  X  X  X 

H^^^  ;=  o  A  0*  =  X  X  X  X  X 

^ToCl  0  Hm\  X  X  X  X 

XXX 
X  X  . 

into  Hessenberg  form  without  changing  6.  Let  CTq  :=  and  oo  :=  Then 

'  8  V  Oq 

fTfl  X  X  X  X  X  X 

0  X  X  X  X  X  X 

G2{qo)^^‘^  =••  =  X  X  X  X  X 

X  X  X  X 

XXX 
X  X  . 

is  a  Hessenberg  matrix  with  a  trailing  principal  (m-t- 1)  X  (m-|- 1)  submatrix,  which  is  both  unitary 
and  of  Hessenberg  form.  On  the  completion  of  the  similarity  transformation  of  we  obtain 

r  6  aj 

<70  X  X  X  X  X  X 

X  X  X  X  X  X 

.^(^>G2(qo)=  (g)  X  X  X  X  X  .  (3.1) 

X  X  X  X 

XXX 
X  X  . 

The  circled  element  in  (3.1)  forms  a  “bulge”,  which  is  to  be  chased  down  along  the  subdiagonal  in 
order  to  obtain  a  matrix  of  Hessenberg  form.  Define  G3(ai)  so  that  G3(qi).^^^^G2(q'o)  =•'  H'^^^  is 
a  Hessenberg  matrix.  Then  has  a  bulge,  which  we  annihilate  by  multiplying  from  the 

left  with  G4(a2).  Proceeding  in  this  manner,  we  ultimately  chase  the  bulge  off  the  bottom  of  the 
matrix,  and  obtain  the  Hessenberg  matrix  .ffi'”"*lG^_i(7m-3),  which  is  unitarily  similar  to  H  ^  . 
The  trailing  principal  (m  -t-  1)  x  (m  -i-  1)  submatrix  of  is  uniUrily  similar  to 

a  unitary  Hessenberg  matrix  with  positive  subdiagonal  elements.  The  latter  matrix  is  the  desired 
Hessenberg  matrix  H'm^i  =  H{^[,Y2, .  ■  ■  ,l'm+i)- 

This  procedure  for  adding  a  weight-abscissa  pair  to  Hm,  if  implemented  by  directly  manipulating  the 
elements  of  the  matrices  H^'‘\  for  1  <  ^'  <  m,  would  require  0{m'^)  arithmetic  operations.  However, 


we  note  that  for  each  k  the  trailing  (m+l)x(m+l)  principal  submatrix  of  is  unitary  and  of  Hes- 
senberg  form,  and,  therefore,  is  unitarily  similar  to  a  unitary  Hessenberg  matrix,  denoted  by 
with  positive  subdiagonal  elements.  Hence,  we  can  carry  out  the  similarity  transformations  by  ma¬ 
nipulating  the  Schur  parameters  of  the  matrices  =  /f (7{, 72>  •  •  •  > •  •  •  ilm)- 

This  gives  rise  to  a  method  that  requires  only  0(m)  arithmetic  operations  in  order  to  add  a  weight- 
abscissa  pair  to  Hm  ■  We  refer  to  this  method  as  the  inverse  unitary  QR  algorithm  because  of  its 
relationship  with  the  unitary  QR  algorithm  presented  in  [Gr2]. 

Inverse  Unitary  QR  Algorithm;  adding  a  weight-abscissa  pair  (u^.  A),  where  jAj  =  1,  >  0. 

/  /  0  ,  O 

<^o  ■■=  V^o  + 

Oo  :=  -vfaQ-,  00  :=  cto/ct^; 
for  k  :=  1,2, . . . ,  m 

^'k  ■—  0k-\‘\J(^l  +  |Qfc-i  +  7fcA*^~-a*;_i|-,- 

:=  0k-iK*^k-i  +  7fcA*“^afc_i)/<r^; 

0k  ■■=  0k-i<^kl<rU 

/92 

l^k  ■—  l^k-l^k  -  ^  (^k-ly 

D 

Thus,  the  inverse  uiiitary  QR  algorithm  can  be  used  to  construct  the  unitary  n  x  n  Hessenberg 
matrix  H  =  //(71, . . . ,7„_i,7„)  from  its  eigenvalues  and  the  first  components  of  its  normalized 
eigenvectors  in  O(n-)  arithmetic  operations. 


4,  An  Inverse  Spectral  Problem 


/  w  ^ 

Let  H  =  H{li, . .  .,7n-i,7n)  and  <0nW  =  -TTvi  =  T~V‘-  q  :=  e’"^  for  some  0  <  r  <  2;r, 


ti’n(A) 


A- A* 


and  consider  the  polynomials 


:=  (1  -  o)Ajrfe(A) -h  a^jt(A),  0  <  k  <  n.  (4.1) 

Proposition  4.1.  The  polynomials  {X7e(A)}5J_q  are  the  monic  Szego  polynomials  corresponding 
with  the  Schur  sequence  {o7fc}]J_j. 

Proof.  The  definition  (4.1)  and  Proposition  2.1  yield 


where  ^*:(A)  :=  (1  -  a)A7f*.(A)  4-  Qt/’*(A).  On  setting  Xk{^)  =  ct^fc(A),  we  obtain 
A'fc('^)1  _  [  ttTfcl  [x*:-i('^)l  [a:o(A)1_[i 

.X-t(A)J  [q7*A  1  J  [xit-iC-^)]  ’  [xo(A)J  [1_  ■ 


An  immediate  consequence  of  Proposition  4.1  is  that  each  zero  of  XnC-^)  has  unit  modulus.  When 
o  =  -1  the  polynomials  {x/t('^)}J?=o  known  as  the  Szego  polynomials  of  the  second  kind 
corresponding  with  {7*’}J5_i. 
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Proposition  4.2.  The  zeros  of  Xn(^)  strictly  interlace  the  zerob  of  on  the 

unit  circle.  Moreover,  nL  i(/u/Afc)  =  a. 

Proof.  Let  A  ;=  for  some  0  <  ^  <  2?:,  and  let  a  =:  e'’’  for  some  0  <  r  <  27r.  Then 


4'n(A):= 


An(A) 


=  E 


Uk- 


A  -  QA;t 


t^-n(A)  '  A  -  A, 


=  E 


Uk 


1  —  AAfc  +  Q:(l  —  AAfc) 


k=l 


=  f^2e'- 


!A- A,|2 

/2  Re(eW^(l- AAfe)) 


Uk 


k=l 


|A-A,P 


Thus,  Im(e  =  0  for  0  <  ^  <  2;r.  Let  A*,  =:  e'** .  We  may  assume  that  0  <  di  <  62  < 

■  ■  ■  <  On  <  27r.  Then 


k=l 


k=i  (1  -  cos(0fc  -  ^))- +  sin-((9/..  -  0) 

cos(r/2)  —  cos(r/2)  cos{6k  —  6)  +  sin(r/2)  sin(0;t  —  0) 
k=i  1  -  cos(0t  -  0) 


=  '^u>k  (cos(r/2)  +  sin(r/2)cot((^*  -  0)/2)) . 


*=i 


From  0  <  r  <  27r  and  Uk  >  0  for  all  k,  it  follows  that 


d  1  " 

=  -sin(r/2)^u;fc/sin\(fi'fc  -  6»)/2)  >  0  for  0  ^  Ok,  I  <  k  <  n. 


k=:l 


Thus,  4>„(e'*)  -*  -00  as  0  \  Ok,  and  $„(e’^)  -♦  oc  as  0  /"  This  shows  that  0  — ►  $„(e’*)  has 
precisely  one  zero  in  ]0fc,^fc+i[.  Consequently,  the  zeros  /ij  of  $n(A)  strictly  interlace  the  Aj  on  the 
unit  circle.  The  second  statement  of  the  Proposition  follows  from  the  fact  that  Xni0)/4’n{0)  =  o.  a 

Let  Xk  ='■  and  p*.  =:  e*"*  for  0  <  <  27r.  Then  we  have  —  Ok)  —  r.  We  may 

assume  that  the  arguments  have  been  ordered  so  that 

0  <  <  Ui  <  Oo  <  ■•■  <  V’n-l  <  On  <  I'n  <  +  27r.  (4.2) 

Proposition  4.3.  With  the  above  notation,  and  the  ordering  (4.2),  the  weights  Wfc  are  given  by 


Wfc  =  - 


Proof,  ^^’e  have 


n  sin(((^j  -  0k)/2) 
J _ i=i _ 

n  -  90/2) 

}=i 
j^k 


(1  -  a)A,a;*  =  lim  (A  -  A,)^4^  =  - 

i=t 

j^k 
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since  \„(A)  and  <,'n(A)  are  monic.  Thus, 


(1  -aW,  = 


n  (f"**  -  0(1-^'^*'’  n  sin((;/j  -  0a.)/2) 

-  =  - 


n  {e'«‘  - )  n  { A  -  * ) 

j=i  i=i 

A- 


.  (T3) 


n  sin{(dj  -  Bk)l2) 
j=i 

J^k 


because 

1  -  e*'^  = -2isin(/3/2),  (-4.4) 

and  X)j=i(A'j  ~  ^j)  —  T"-  The  formula  for  the  weights  now  follows  by  substituting  (4.4)  with  P  :=  r 
into  (4.3).  □ 

We  can  now  state  the  inverse  spectral  problem  and  its  solution. 

Problem:  Given  two  sets  of  n  mutually  interlacing  points  on  the  unit  circle  {A*..}^_j  and 

determine  the  unique  unitary  llessenberg  matrix  H  =  . 'yn-i<'>n)  ^^nd  |a|  =  1  such  that 

the  spectrum  of  H  is  ■(At}’J_,  and  the  spectrum  of  . . . ,  Q-)n_i ,  a*;n )  is 

Solution:  Let  a  ;=  n”=i(/'A/A/c),  calculate  the  weights  by  Proposition  4.3,  and  use  the  inverse 
unitary  QR  algorithm  to  construct  II  =  II{~i\ . 'Jn-iOn)- 
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